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3? (57) Abstract: The invention concerns a methodfor separating microseismic signals from seismic signals acquired in the context of 
2! active seismic monitoring operations of underground zones being exploited. It consists in producing seismic recordings from signals 
^5 emitted by a single or several seismic sources monitored by signals orthogonal relative to one another. In that case the signals are 
processed to separate the respective contributions of the seismic sources to the received signals and reconstruct the seismograms 
equivalent to those which would be obtained by separately actuating the seismic sources. The induced microseismic signals (passive 
seismics) are separated from the seismic signals derived from sources 5 (active seismics) essentially by isolating the contribution of 
fs| the latter, by comparison with a reference spectral model to the fundamental frequencies emitted and to their respective harmonics, 
and by reconstructing in the time domain the induced microseismic signals and the signals derived from the seismic sources 5. The 



invention is applicable to seismic monitoring of underground hydrocarbon or fluid storage reservoirs. 
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(57) Abrege : MSthode pour separer des signaux de microsismicite' induite de signaux sismiques acquis dans le cadre d'operations 
de surveillance sismique active d'une zones souterraines en cours d' exploitation. On forme des enregistrements sismiques a partir 
de signaux emis par une source sismique unique ou plusieurs sources sismiques pilotees par des signaux orthogonaux les uns rela- 
tivement aux autres. Dans cette alternative les signaux sont traites pour separer les contributions respectives des sources sismiques 
aux signaux recus et reconstruire les sismogrammes equivalents a ceux que Ton obtiendrait en actionnant separement les sources 
sismiques. On separe les signaux de microsismicite' induite (sismique passive) des signaux sismiques 6manant des sources 5 (sis- 
mique active) essentiellement en isolant la contribution de ces derniers, par comparaison avec un modele spectral de reference aux 
frequences fondamentales 6mises et a leurs harmoniques respectives, et en reconstituant dans le domaine des temps les signaux de 
microsismicite^ induite et les signaux Smanant des sources sismiques 5. - Applications a la surveillance sismique de reservoirs sou- 
terrains d'hydrocarbures ou de stockage de fluides. 
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METHODE DE SEPARATION DE SIGNAUX DE MICROSISMICITE DE 
5 SIGNAUX SISMIQUES EMIS PAR UNE OU PLUSIEURS SOURCES 



L'invention concerne une m6thode pour discriminer des signaux microsismiques 
induits par Texploitation d'un gisement souterrain, parmi des signaux sismiques 6mis dans 
le cadre d'operations de surveillance sismique 

10 Pour suivre Involution d'un gisement souterrain en cours Sexploitation, il est connu 

d'implanter k poste fixe des sources sismiques et des r^cepteurs sismiques, en surface ou 
dans des puits, qu'il s'agisse de puits d'exploration ou Sexploitation du gisement et de 
proc6der k intervalles de temps fixes k des operations d'enregistrement des signaux 
sismiques 6manant de la formation. II peut s'agir de signaux sismiques induits par 1'activite 

15 d'exploitation du site ou de signaux sismiques renvoy6s par les discontinues geologiques 
de la formation, en r^ponse k remission d'6branlements dans le sol au moyen d'une ou 
plusieurs sources sismiques. La comparaison des enregistrements effectu6s k des instants 
diffcrents de Texploitation du site, fournit des donn^es sur son evolution. 

Differents systemes de surveillance sismique k long terme sont d6crits par exemple 
20 dans les brevets EP 591 037 (US 5 461 594), FR 2 593 292 (US 4 775 009), FR 2 728 973 
(US 5 724 311) ou FR 2 775 349. 

Par le brevet EP 748.457 (US 5 724 311) de deux des demandeurs, on connait un 
systeme de surveillance permanente permettant d ! obtenir une parfaite reproductibilit6 des 
conditions op6ratoires dans des operations r6petitives de surveillance sismique d'une zone 
25 souterraine travers6e par au moins un puits ou forage, et notamment d'un reservoir 
souterrain de stockage de gaz. Ce systeme comporte, installs k poste fixe, un ou plusieurs 
ensemble(s) de r6cepteurs sismiques (enterr6s en surface ou places dans un ou plusieurs 
puits), plusieurs sources sismiques r6p6titives (enterrees ou en surface) ainsi que d ! un 
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rtseau de liaison permanente pour Talimentation selective de ces sources en tnergie. Une 
station centrale teltcommande stlectivement chacune de ces sources, et enregistre les 
signaux sismiques tmanant de la zone souterraine en rtponse aux ondes sismiques 
transmises stlectivement dans le sol par les sources. 

5 GrSce k cet ensemble de sources k poste fixe dont le couplage avec les terrains 

environnant reste stable, et k ce rtseau d'alimentation au moins en partie enterrt et dont 
1'aire d'emprise en surface est rtduite, on peut mener toute une strie d'optrations sismiques 
de surveillance de longue duree dans des conditions operatoires stables, sans risque 
^incompatibility avec les activites du chantier Sexploitation. 

10 Par le brevet FR 2.775.349 (US 6 182 012) des deux memes demandeurs, on 

connait une mtthode et un dispositif de surveillance permanente d ! une formation 
souterraine dans laquelle on forme des cavitts pour des sources sismiques (de preference 
assez profondes pour qu'elles soient coupltes avec la formation au-dessous de la zone 
alttrte (WZ)), on fore au moins un puits soit au-dessous de chacune des sources ou au 

15 voisinage immediat de chacune d'elles. Dans chacun des puits ainsi forts, on descend une 
antenne constitute d'une pluralite de recepteurs sismiques (geophones, hydrophone, etc.) 
que Ton couple avec la formation environnante. Avec ce dispositif, on rtalise des cycles 
d'tmission-rtception avec Emission Sondes sismiques dans la formation par au moins une 
source et acquisition des signaux captts par les recepteurs sismiques, en rtponse aux ondes 

20 renvoy6es par la formation. Les recepteurs peuvent par exemple etre associ6s avec un tube 
de transfert de fluides descendu dans le puits, dans le but de relier une zone souterraine k 
un appareil Sexploitation en surface et coupl6s avec la formation environnant le puits. 

Les differentes sources du dispositif de surveillance sismique peuvent etre 
actionn6es successivement, en m6nageant entre les declenchements un intervalle de temps 
25 suffisant pour recevoir les ondes renvoyees par la zone investigu6e. On peut aussi utiliser 
plusieurs sources sismiques tmettant les memes signaux que Ton declenche simultan6ment 
pour augmenter la puissance 6mise. 

Par le brevet FR 2.805.051 des demandeurs, on connait une m6thode et un syst&me 
de surveillance sismique d'une formation souterraine. Comme schematises aux Fig.l k 4, le 
30 systeme comporte par exemple un r6seau d'antennes sismiques 2 constitutes chacune d'un 
ensemble de capteurs sismiques 4 disposes k intervalles reguliers le long d'un puits 3 fort 
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dans le sol. Ce r£seau peut §tre rdgulier comme sch6matis6 sur la Fig.2, ou irr6gulier. Les 
capteurs sont par exemple des g6ophones mono-directionnels orientes verticalement ou 
multi-axes (triphones) et/ou des hydrophones. A proximity de chaque antenne 2, est 
disposee une source sismique 5. Comme sources, on utilise avantageusement des vibrateurs 

5 de type pi6zo61ectrique, tels que d6crit dans le brevet FR 2.791.780 (US 6.338.394) aux 
noms conjoints des demandeurs, qui sont installs k demeure au voisinage imm6diat de 
chaque antenne 2. A chaque antenne peut Stre associe un module electronique de contrdle 
et d'enregistrement 6 (Fig.l). Ces modules electroniques 6 peuvent Stre reli6s k un 
laboratoire central de commande et de synchronisation 8. Suivant un autre mode de mise en 

10 ceuvre, sch6matis6 sur la figure 4, les antennes 2 sont relives directement k un laboratoire 
central 8 regroupant tous les modules 61ectroniques individuels de contr61e et 
d'enregistrement 6 (Fig.l). 

Des ondes sismiques sont 6mises dans la formation souterraine par une ou plusieurs 
sources sismiques (vibrateurs) 5. Dans le cas de sources multiples, celles-ci 6mettent 

15 simultan6ment ; elles sont pilot6es par des signaux orthogonaux les uns relativement aux 
autres, de maniere k former un signal vibratoire composite. Les vibrateurs fonctionnent en 
permanence ou tout du moins sur une dur6e assez longue pour obtenir un rapport signal au 
bruit suffisant. Les ondes sismiques gen6r6es par les sources sismiques 5 se propagent vers 
le bas (ondes descendantes 9). Ces ondes incidentes sont tout d'abord enregistr6es par les 

20 r^cepteurs 4 de chaque puits 3. Les ondes renvoy6es par les discontinues de la zone 
(interfaces sismiques) se propagent vers le haut. Ces ondes montantes 10 sont enregistr6es 
aussi par les differents r6cepteurs 4. De la sorte, les ondes montantes et descendantes sont 
superpos6es sur les sismogrammes. On applique aux enregistrements les traitements usuels 
bien connus des gens de Tart convenant k ceux obtenus par. la. technique dite de PSV 

25 (Profils Sismiques Verticaux). Ensuite, on discrimine les contributions respectives des 
sources sismiques au signal vibratoire composite et on reconstruit les sismogrammes 
Equivalents k ceux que Ton obtiendrait en actionnant s6parement les sources sismiques. 
Comme signaux orthogonaux, on utilise par exemple des signaux sinusoidaux de 
frequences differentes les unes des autres aussi bien par leurs composantes fondamentales 

30 que par leurs harmoniques respectives, ou des signaux form6s k base d'ondelettes, de 
polynomes de Legendre ou de series al6atoires, etc. 
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Dans le cas notamment oil les signaux orthogonaux 6mis sont des sinusoi'des, on 
realise par exemple la discrimination des contributions respectives des sources sismiques, 
par determination de T amplitude et de la phase du signal vibratoire composite aux 
frequences fondamentales des signaux de pilotage appliqu6s aux sources sismiques. Cette 

5 discrimination comporte par exemple une pond6ration des signaux enregistr€s par un 
facteur de pond6ration (ou d'apodisation) en cloche et une determination de 1' amplitude et 
de la phase du signal composite, ou bien encore une selection par transformee de Fourier, 
de raies du spectre complexe associees respectivement aux differents signaux ponder6s. La 
reconstruction des sismogrammes correspondant specifiquement aux differentes sources 

10 sismiques est effectu6e par exemple en appliquant, apres leur separation, une transformee 
de Fourier inverse, aux raies associees respectivement aux differents signaux pond6res. 

Par ailleurs, il est connu que l'exploitation des reservoirs induit des variations de 
contraintes mecaniques dans le reservoir et dans sa couverture. Ces variations de contrainte 
ont pour effet de generer des signaux sismiques se propageant dans les terrains. Les 
15 recepteurs des differentes antennes 2 vont enregistrer ces signaux (ondes montantes 11). 

Comme les sources sismiques 5 sont situees beaucoup plus pr£s des recepteurs 4 
des antennes 2 que du reservoir et que renergie qu'elles emettent est bien plus intense que 
renergie des signaux microsismiques induits generes dans le reservoir, il en resulte que les 
signaux microsismiques sont generalement masques et indecelables. 

20 La methode selon l'invention 

La m6thode selon l'invention permet de separer des signaux de microsismicite 
induite et des signaux sismiques multi-sources dans le cadre d'operations de la surveillance 
d ! une zone souterraine en cours d'exploitation. 

Elle comporte remission d* ondes sismiques dans la formation en couplant avec elle 
25 une ou plusieurs sources sismiques. Dans ce dernier cas, celles-ci emettent simultanement 
des signaux orthogonaux les uns relativement aux autres, de manure k former un signal 
vibratoire composite, la reception des signaux renvoyes par la formation en reponse k 
remission des ondes sismiques, Tenregistrement des signaux re9us par au moins un 
capteur sismique et le traitement des signaux enregistres pour separer les contributions 
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respectives des sources sismiques aux signaux re§us et reconstruire les sismogrammes 
Equivalents k ceux que Ton obtiendrait en actionnant s6parement les sources sismiques. 

On separe les signaux de microsismicite induite (sismique passive) des signaux 
sismiques 6mis essentiellement en isolant la contribution de chacune des sources par 
5 comparaison avec un module spectral de reference aux frequences fondamentales emises et 
k leurs harmoniques respectives, et en reconstituant dans le domaine des temps les signaux 
de microsismicite induite et les signaux sismiques 6mis. 

La contribution spectrale des signaux de microsismicite au spectre des signaux 
re§us est obtenue par exemple en retranchant les valeurs Samplitude et de phase associees 
10 au modele spectral de reference, de T amplitude et de la phase du spectre associe aux 
enregistrements. 

Suivant un mode de mise en ceuvre, comme modele spectral de reference, on utilise 
par exemple un module courant forme par mise k jour d'un modele spectral prenant en 
compte la contribution de cycles d'enregistrement precedents. 

15 Le module spectral courant peut etre forme, en determinant une valeur moyenne des 

spectres de frequence formes k partir d'enregistrements (anterieurs et/ou posterieurs) 
obtenus pour la mSme source et les memes frequences. 

Suivant un autre mode de mise en oeuvre, on forme le modele spectral courant, en 
determinant une valeur mediane des spectres de frequence formes k partir d'enregistrements 
20 anterieurs obtenus pour la meme source et les m§mes frequences. 

Suivant un autre mode de mise en oeuvre, on forme le modele spectral courant, par 
extrapolation ou interpolation k partir du spectre de frequence k partir de valeurs spectrales 
voisines. 

Suivant un autre mode de mise en oeuvre, on forme le modele spectral courant en 
25 utilisant la raie obtenue k cette frequence. 

Suivant un mode de mise en ceuvre, la discrimination de signaux de microsismicite 
induite parmi des signaux emis dans le cadre d'operations de surveillance sismique active 
est obtenue par les 6tapes suivantes : 
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a) pour chaque enregistrement n d\in cycle p d'enregistrement, on calcule les 
contributions respectives des differentes sources aux frequences fondamentales ; 

b) on calcule ensuite le rapport de la contribution k un modele spectral courant 
forme par actualisation d'un modele spectral precedent k partir de frequences 6mises durant 

5 Tenregistrement precedent ainsi que de leurs harmoniques ; 

c) on deduit la partie de V enregistrement n du cycle p, que Ton peut associer aux 
op6rations de surveillance sismique active ; 

d) on deduit la partie de F enregistrement n du cycle p, relevant proprement de 
l'activite microsismique passive ; et 

10 e) on forme par inversion dans le domaine des temps, des contributions spectrales 

respectives de chacun des sources sismiques aux frequences fondamentales et k leurs 
harmoniques, apres completion d'un cycle de mesures, les sismogrammes que Ton peut 
associer aux operations de surveillance sismique active ; et 

f) on forme par inversion dans le domaine des temps, k partir de la partie relevant 
15 de l'activite microsismique passive, les signaux microsismiques sous-jacents contenus dans 
les enregistrements. 

La contribution spectrale est obtenue en multipliant une fonction de transfert entre 
une ondelette caracteristique de la source et un sismogramme associe a un recepteur donne, 
par une ondelette caracteristique de la source. Cette fonction de transfert peut etre 
20 continuellement mise k jour au cours d ! un cycle courant k partir d'une estimation faite au 
cours d ! un cycle precedent et d f une estimation initiale faite au cours du cycle courant par la 
relation. 

Presentation sommaire des figures 

Les caracteristiques et avantages de la methode selon Tinvention, apparaitront plus 
25 clairement k la lecture de la description ci-apres d'un exemple non limitatif de mise en 
oeuvre, en se referant aux dessins annexes oii : 



la Fig. 1 montre schematiquement le dispositif de surveillance mis en place sur le terrain 
pour surveiller sismiquement une formation souterraine ; 
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- la Fig.2 montre un mode de disposition des unites demission-reception sur le teirain ; 

- la Fig.3 montre schematiquement une de ces unites d ! 6mission-reception comprenant 
une source sismique et une antenne de reception ; 

- la Fig.4 montre une variante du dispositif de surveillance de la Fig. 1 ; et 

5 - la Fig.5 montre des exemples de signaux et spectres de frequence mentionnes dans le 
cours de la description. 

Description detaill£e 

La m&hode permet done de r^aliser des operations de surveillance sismique d'une 
zone souterraine en utilisant un ensemble de capteurs sismiques et une ou plusieurs sources 
10 sismiques vibratoires. Dans le cas d'utilisation de plusieurs sources simultanees, celles-ci 
sont actionn6es simultan6ment par des signaux k des frequences diff6rentes choisies de 
fa?on que Ton puisse discriminer les contributions de chaque source sur les sismogrammes 
obtenus & partir des signaux requs et enregistres. 

Ceci est realise d'une fagon generate par le pilotage des diff6rentes sources par des 
15 sinusoides de frequences differentes et par utilisation de techniques de calcul num6rique 
connues telles que la transformee de Fourier inverse pour sSparer les contributions aux 
sismogrammes obtenus des differentes sources sismiques. 

La distorsion ne peut 6tre negligee avec des sources m£caniques. En raeme temps 
qu'une frequence ft, une meme source Si emet les frequences 2fi, 3ft ... nft. En 
20 consequence, si ft et fj sont les frequences respectives des deux sources Si et Sj de 
l'ensemble de sources, il convient si on veut separer leurs contributions respectives non 
seulement que £ ^ fj mais aussi que ft * 2fj, ft * 3fj , . . . ft * nfj . 

Comme d6jk decrit dans le brevet FR 2.805.051 precite, si on emet le signal 
composite P t constitue de la somme de N sinusoides {ft, Ai, <I>, } avec 1 < i < N, toutes les 

25 frequences etant contenues dans une bande spectrale comprise entre deux frequences 
limites ft, et ff, le sismogramme T t observe au point R aura pour transformee de Fourier & la 
frequence ft,, le nombre de module Ai et de phase <&, 6gaux & 1' amplitude et k la phase de 
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la sinusoi'de TV On peut ainsi, en emettant successivement toutes les sinusoides de 
frequences fb k ff, reconstruire par transform6e de Fourier inverse le sismogramme T t . 

Dans le cas oil, par exemple, toutes les amplitudes Ai sont egales k 1 et toutes les 
phases sont nulles, le signal P t obtenu est tres proche de celui resultant de Pauto- 

5 conflation d'un signal k frequence glissante variant dans Pintervalle [fb - ff,] (sweep), 
utilise couramment en vibrosismique. D'aprSs la th6orie de la transform6e de Fourier 
discrete, bien connue des gens de Tart, si Ton desire ecouter la source Si pendant le temps 
te, P increment de frequence entre les sinusoides est Af = 1/te et le nombre de sinusoides 
n6cessaires est Nf = (ff - f b )te. 

10 On peut done exciter simultan&nent N vibrateurs SI, S2, Si,... Sn installs sur le 

terrain au moyen de signaux vibratoires de frequences telles que chaque source est excise 
successivement par chacune des Nf sinusoides ci-dessus k chaque instant, sous reserve que 
les frequences respectives des sinusoides emises k un meme instant par les diff6rentes 
sources sismiques sont toutes differentes les unes des autres. On realise ainsi avec chaque 

15 vibrateur Nf enregistrements correspondant aux Nf frequences seiectionnees dans 
Tintervalle de balayage. La separation des signaux re?us par les capteurs sur le terrain, en 
reponse k P emission simultanee des differents signaux est ainsi obtenue par selection de la 
raie k la frequence appropriee. 

On repete ainsi des cycles de Nf enregistrements. L'enregistrement En d'un cycle 
20 d' acquisition donne contient un ensemble de frequences fondamentales et harmoniques qui 
ne depend que de n. Partant de cette observation, la methode que Ton va appliquer pour 
detecter les effets de la microsismicite, repose sur le constat verifie dans la pratique, que les 
variations d' amplitude et de phase observees k la frequence f et k Pinstant t sont trfes peu 
differentes de celles observees k la frequence / — Af et a Pinstant t-At. Cela permet de 
25 former un modele spectral cumulatif construit k partir d'enregistrements de cycles 
precedents aux memes frequences mais forme done sur un temps long pendant lequel les 
caracteristiques des signaux sismiques emis sont susceptibles de varier, et adapte en 
permanence k partir d'enregistrements obtenus pendant une periode suffisamment courte 
pour que les variations de signal puissent etre negligees L'efficacite de la methode 
30 proposee est optimale quand les frequences successives emises par une m6me source 
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s'incrementent r6gulierement. On pourra, par exemple, alterner des cycles en frequences 
montantes et en frequences descendantes. 

Pour parvenir dans la pratique & la mise en evidence de cette contribution spectrale 
des signaux microsismiques sous-jacents, on va proc6der comme suit. 

5 Notations : 

■ fi,n>va designer la frequence fondamentale emise par la source i pendant 
l'enregistrement n (cette frequence est independante du cycle) ; 

Pour simplifier les notations des grandeurs suivantes, Tindice designant le numero 
du recepteur est omis. Quel que soit ceku-ci, les grandeurs sont calcuiees de la meme 
10 manure. 

■ R Pt n va designer l'enregistrement de rang n du cycle p ; 

■ Cp^u, la contribution de la source i a renregistrement n du cycle p ; 

■ Mp, n , le- modele initial de la partie active de l'enregistrement n du cycle p. Ce 
module ne contient que les frequences emises ainsi que leurs haimoniques pendant 

15 renregistrement n ; 

■ Ep >n , le rapport entre les contributions actives de l'enregistrement n du cycle p et 
le modele initial de la partie active de l'enregistrement n du cycle p ; 

■ A p>n , la partie active de l'enregistrement n du cycle p ; 

■ P p>n , la partie passive de l'enregistrement n du cycle p (microsismique) ; 

20 ■ k le coefficient d'actualisation des ecarts (typiquement de l'ordre de 15 k 

25%) ; 

■ h le coefficient d'actualisation des modeles (typiquement de l'ordre de 5 k 10%). 

Compte tenu de ces notations, la methode consiste dans un premier temps, pour 
chaque enregistrement n du cycle p, au calcul des contributions respectives Cp,i >n des 
25 differentes sources aux frequences fondamentales ainsi qu'k leurs harmoniques par la 
methode decrite dans le brevet FR 2.805.051 precite. Du fait de la presence eventuelle 
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d'energie microsismique, la contribution des diff&rentes sources est susceptible d'6tre 
bruit6e. 

On peut am61iorer cette estimation en associant les contributions de tous les 
r6cepteurs d'un meme enregistrement C p ,i, n . Ces contributions peuvent etre consid6r6es 
5 comme le produit de deux facteurs dont Tun ne d6pend que de la source (et non du 
r6cepteur) et Tautre du r6cepteur et non de la source. 

Si l'indice r repr6sente le r6cepteur, 
C = S T 

T ne depend de p, i et n que par la variation des caracteristiques du rScepteur et de 
10 son environnement avec le temps. Ces variations sont trds lentes et peuvent etre negligees 
ou estim6es de fa§on plus precise. 

S et T sont d6termin6s k un facteur pres. On peut ainsi choisir pour la moyenne 
des ondes descendantes. Cette moyenne est obtenue par un calcul bien connu de ceux qui 
traitent les PSV (mise k plat des ondes descendantes et moyenne). Les T p ,i, n , r sont alors les 
15 fonctions de transfert entre Fonde descendante et le sismogramme. 

Dans le cas oil Ton d6sire prendre en compte la variation de T, on peut utiliser par 
exemple le coefficient d'actualisation h pour mettre k jour k partir de T PfitI1 .i tr en 

procedant de la mani&re suivante : 

Si Co p ,i, n ,r d6signe F estimation initiale de la contribution C et To Pt i, n ^ celle du 
20 facteur r6cepteur, on calcule S par mise k plat et moyenne des ondes descendantes puis 
To P) i >n , r par la relation : 

To • 

A ^p,i,n,r q 

On obtient F estimation finale du facteur r6cepteur par : 

T p ,i, n ,r = (l~ h ) T p.U-l,r + hTO p4>n>r 



25 



et la contribution finale de la source i au cycle p de F enregistrement n par : 
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P — Q T 

^p.i.n.r ^p,U x p,i.n,r 

On calcule ensuite le rapport Ep^ de P enregistrement courant au modele spectral 
par la relation : 

1V1 P,n 

5 Ce calcul n6cessite la connaissance de M p>n dont le calcul d6crit ci-dessous est 

effectu6 lors du cycle p-L D comporte un lissage sur les frequences prec6dentes de fa?on a 
limiter Peffet de P6nergie microsismique dans son estimation. La formule propos6e 
n'utilise que les enregistrements pr£c6dents et autorise done le temps r6el. Dans les cas 
difficiles, on peut remplacer cette somme par une moyenne ou une m6diane courante. II 

10 faut alors travailler en temps differe (le temps n6cessaire h r6aliser la seconde moitie des 
enregistrements utilises dans le calcul de la moyenne). 

Si M p+U repr6sente le modele initial de la partie active de P enregistrement n du 
cycle p+1, alors : 

M p ^=a-h)M Ptn +li^E p<n _ 1 C p ^ n _ 1 

15 Ce module consiste essentiellement en une moyenne ponder6e de la contribution 

courante et des contributions pr6c6dentes. La formule utilis6e autorise la variation du 
modele : la contribution du cycle p-q au module utilise pour le cycle p est multipliee par le 
coefficient (l-h) q qui tend vers z6ro d'autant plus vite que h est voisin de 1. Cela permet au 
modele de suivre les variations saisonnieres, et aux rapports E p ^ de ne refl&er que les 

20 variations journalieres.c) Connaissant le modele spectral Mp >n ainsi que le rapport Ep, n de la 
contribution au modele pr6c6demment etablie, on en ddduit la partie A PtI1 de 
P enregistrement n du cycle p, relevant de la sismique active Partie active : 

De la meme fagon la partie P p , n de P enregistrement n du cycle p, relevant 
25 proprement de Tactivit6 microsismique (passive) se calcule par la relation : 



Pp,n — ^-p,n — -^p»n 
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Les contributions spectrales respectives Cp,i, n de chacune des sources Si aux 
frequences fondamentales et & leurs harmoniques, apr&s completion d'un cycle de mesures 
fournissent par FFT inverse, les sismogrammes (dans le domaine des temps) utilises dans 
le cadre de r exploitation active. 

5 Par FFT inverse, h. partir de la partie P p , n relevant de contributions passives, on 

obtient V estimation de l'activite microsismique sous-jacente contenue dans 
Tenregistrement R p , n . 

La mise & jour du rapport Ep >n entre la contribution courante et le modele spectral 
M p , n , est effectuee de fagon relativement rapide de fagon a permettre des variations rapides 
10 du signal. La mise k jour du modele initial s'effectue de fagon beaucoup plus lente ; en fait 
le modele initial pourrait ne pas etre mis k jour. 

L'interet de ces mises & jour est d'optimiser Testimation du signal actif quand ses 
variations ne peuvent etre controlees (variations saisonnieres et journalieres). Dans les cas 
les plus favorables, elles ne sont pas n^cessaires et Testimation de la partie active se 
15 r6duira au modele initial. 

Une autre simplification possible est de n6gliger le bruit sismique aux frequences 
6mises (et, si besoin, aux harmoniques de ces frequences). On ne se preoccupe plus 
d'estimer la partie active et on se contente, pour obtenir la partie passive, d'annuler ou 
d'interpoler les raies correspondant aux frequences emises (et, s'il le faut, celles de leurs 
20 harmoniques). 

Dans un premier temps, on supprime les raies harmoniques du spectre de frequence 
de chaque signal enregistre et on interpoleT amplitude et la phase aux valelirs de raie & la 
frequence de la source ; et puis on reconstitue par inversion Tenregistrement non affecte 
par les emissions des sources contenant eventuellement un signal microsismique. Ce mode 
25 de mise en oeuvre convient plutot dans le cas d'un nombre de sources peu eieve et/ou 
pr6sentant une variabilite importante du signal avec le temps. 
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REVENDICATIONS 



1) M6thode de surveillance sismique active d'une formation souterraine (1) 
permettant une separation de signaux de microsismicit^ induite de signaux sismiques 6mis 

5 dans le cadre d'operations de surveillance sismique active d'une zone souterraine en cours 
Sexploitation, comportant la realisation de cycles d'enregistrement sismique avec emission 
d'ondes sismiques dans la formation en couplant avec elle une ou plusieurs sources 
sismiques (5) qui, dans ce cas, 6mettent simultan6ment des signaux orthogonaux les uns 
relativement aux autres, de mani&re k former un signal vibratoire composite, reception des 

10 signaux renvoy£s par la formation en reponse k remission des ondes sismiques, 
enregistrement des signaux re?us par au moins un capteur sismique (2) et traitement des 
signaux enregistr^s pour s^parer les contributions respectives des sources sismiques aux 
signaux regus et reconstruire les sismogrammes Equivalents a ceux que Ton obtiendrait en 
actionnant separ6ment les sources sismiques, caract€ris6e en ce que Ton s6pare sur les 

15 enregistrements les signaux de microsismicite induite des signaux sismiques resultant des 
operations de surveillance active, en isolant leur contribution par comparaison avec un 
module spectral de reference en tenant compte des contributions spectrales de chaque 
source (Si) aux frequences fondamentales emises et a leurs harmoniques respectives, et en 
reconstituant par inversion dans le domaine des temps les signaux de microsismicite. 

20 2) Methode selon la revendication 1, caract6ris£e en ce que la contribution spectrale 

des signaux de microsismicite au spectre des signaux regus est obtenue en retranchant les 
valeurs Samplitude et de phase associ6es au modele spectral de reference, de F amplitude et 
de la phase du spectre associe aux enregistrements. 

3) Methode selon la revendication 1 ou 2, caract6ris6e en ce que le modele spectral 
25 de reference est un modfele courant forme par mise k jour d ! un modele spectral precedent en 

tenant compte de la contribution spectrale de cycles d f enregistrement precedents. 

4) Methode selon Tune des revendications pr6cedentes, caracterisee en ce que Ton 
forme le module spectral courant en determinant une valeur moyenne des spectres de 
frequence formes k partir d ! enregistrements anterieurs et/ou posterieurs obtenus pour la 

30 meme source et les memes frequences. 
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5) Methode selon Tune des revendications pr6c6dentes, caract6ris6e en ce que Ton 
forme le moddle spectral courant en determinant une valeur m6diane des spectres de 
frequence formes h partir d'enregistrements anterieurs et/ou posterieurs obtenus pour la 
meme source et les memes frequences. 

5 6) M6thode selon Tune des revendications pr6c6dentes, caract6ris6e en ce que Ton 

forme le module spectral courant par extrapolation ou interpolation & partir du spectre de 
frequence & partir de valeurs spectrales voisines. 

7) Methode de surveillance sismique active d'une formation souterraine permettant 
une discrimination de signaux de microsismicite induite parmi des signaux emis dans le 

10 cadre d'opSrations de surveillance sismique active d'une zone souterraine en cours 
Sexploitation, comportant la realisation de cycles d f enregistrement sismique avec emission 
d'ondes sismiques dans la formation en couplant avec elle N sources sismiques (Si) 
emettant simultan^ment et pilotes par des signaux orthogonaux les uns relativement aux 
autres, de mani&re k former un signal vibratoire composite, reception des signaux renvoy6s 

15 par la formation en r6ponse & remission des ondes sismiques, enregistrement des signaux 
re$us par des moyens de reception sismique et traitement des signaux enregistres pour 
sdparer les contributions respectives des sources sismiques aux signaux regus et 
reconstruire les sismogrammes Equivalents k ceux que Ton obtiendrait en actionnant 
s6par6ment les sources sismiques, caract6ris6e en ce qu'elle comporte les etapes suivantes : 

20 a) pour chaque enregistrement n d'un cycle p d'enregistrement, on calcule les 

contributions respectives (C Pt i, n ) des differentes sources aux frequences fondamentales ; 

b) on calcule ensuite le rapport (Ep, n ) de la contribution k un module spectral (M p , n ) 
courant forme par actualisation d f un modele spectral precedent (M p , n .i) k partir de 
frequences emises durant 1'enregistrement (n-1) precedent ainsi que de leurs harmoniques ; 

25 c) on deduit la partie (A p>n ) de T enregistrement n du cycle p, que Ton peut associer 

aux operations de surveillance sismique active ; 

d) on deduit la partie (P p ^) de 1' enregistrement n du cycle p, relevant proprement de 
Tactivite microsismique passive ; et 
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e) on foime par inversion dans le domaine des temps, des contributions spectrales 
respectives (C^ de chacun des sources sismiques (Si) aux fr6quences fondamentales et k 
leurs harmoniques, aprfes completion d'un cycle de mesures, les sismogrammes que Ton 
peut associer aux operations de surveillance sismique active ; et 

5 f) on fonne par inversion dans le domaine des temps, k partir de la partie (P p , n ) 

relevant de Tactivit6 microsismique passive, les signaux microsismiques sous-jacents 
contenus dans les enregistrements. 

8) Methode selon la revendication 7, caracteris^e en ce que la contribution spectrale 
(C Pfi , n ) est obtenue en multipliant une fonction de transfert (T Pfi , n , r ) entre une ondelette 

10 caract6ristique de la source et un sismogramme associe au r^cepteur r, par une ondelette 
caract6ristique de la source. 

9) M6thode selon la revendication 8, caract&isee en ce que Ton met 
continuellement k jour la dite fonction de transfert. 

10) M6thode selon la revendication 9, caracterisee en ce que la mise k jour de la dite 
15 fonction de transfert (T p>i , n ^ ) est obtenue au cours d f un cycle courant k partir tfune 

estimation (Tp^i^) faite au cours d ! un cycle pr6cedent et d f une estimation initiale (To p>ifI1>r 
) faite au cours du cycle courant par la relation : 

W=( 1 - h ) T p,U.U+hT0 p4(ntr 




FIG.4 



WO 2004/034090 



7FR2003/002874 



3/3 



FIG.5 



p,n 



R P,n 


i 




C P,I,n 






1 1 i 1 i,i 




i 




M P,n 


1 




1 1 . 


■ 








Mil 1 1 III 1 1 1 


Ap,n 




1 


1 1 i i . . ■ 


1 — 







INTERN 



NAL SEARCH REPORT 



Inter n^^vllcatlon No 

PCl7h^R/02874 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 G01V1/28 



According to International Patent Classification ((PC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 G01V 



Documentation searched other than minimum documentation to the extent that such documents are included In the fields searched 



Electronic data base consulted during the International search (name of data base and, where practical, search terms used) 

EPO-Internal , WPI Data, INSPEC, COMPENDEX 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category c 



Citation of document, with Indication, where appropriate, of the relevant passages 



Relevant to claim No. 



US 6 442 489 Bl (SOKOLIN HAIM 
27 August 2002 (2002-08-27) 



ET AL) 



column 3, 
column 4, 



line 31 - line 48 
line 40 - line 54 



column 5, line 37 -column 6, line 25 

FR 2 805 051 A (6E0PHYSIQUE CIE GLE) 
17 August 2001 (2001-08-17) 
cited 1n the application 
claim 1 



US 5 377 104 A (S0RRELLS GORDON G 
27 December 1994 (1994-12-27) 
abstract 



ET AL) 



US 6 049 508 A (DEFLANDRE JEAN-PIERRE) 
11 April 2000 (2000-04-11) 
claim 1 



1,7 
1,7 



□ 



Further documents are listed in the continuation of box C. 



ID 



Patent family members are listed In annex. 



° Special categories of cited documents : 

•A" document defining the general state of the art which Is not 

considered to be of particular relevance 
"E* earlier document but published on or after the international 

filing date 

"L" document which may throw doubts on priority claim(s)or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

a O* document referring to an oral disclosure, use, exhibition or 
other means 

*P a document published prior to the International filing date but 
later than the priority date claimed 



T later document published after the International filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

•X' document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

•V document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
In the art 

document member of the same patent family 



Date of the actual completion of the International search 



19 February 2004 



Date of mailing of the International search report 



26/02/2004 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentiaan 2 
NL-2280 HV Rljswijk 
Tel. (+31-70) 340-2040. Tx. 31 651 epo nl, 
Fax: (431-70) 340-3016 



Authorized officer 



Lome, B 



Form PCT/1SA/210 (second shoot) (July 1992) 



INTERNATIONAL SEARCH REPORT 



Inter^ 

PCi/ 



ipllcation No 

1/02874 



Patent document 




Publication 




Patent family 




Publication 


cited in search report 




date 




member(s) 




date 


HQ AAAOA&Q 


R1 


97— nft— poop 


IV U IMC 








r i\ touDuoji 


A 
ft 


1 7— nft— 9001 


CD 

r i\ 


pfto^om 


A1 

MX 


x/ uo tuui 








RR 

Dl\ 


U1UHHOO 


A 
ft 


oa-oi -?oop 

UO UX LuUL 








r a 

Uft 




A1 

ft X 


Ifi-Oft-POOl 

1U UO LUUl 








rw 

v^lM 


X ODOUHD 


T 
i 


07-0ft-P00P 

U/ UO tUvt 








C7 


onni vzop 


A^ 

ftO 


1 ?-0P-P00P 








FA 
Lft 




R1 

DX 


?fi_1 9-900P 








PP 


1 1 7**781 
XI/ »>/ ox 


Al 
ftX 


90_ni -POOP 
£0 UX lUUl 








wu 


01 RQdftl 


A1 
ft X 


Ifi-Oft-POOl 

1U UO £UU X 








HII 


0P0PRQ1 


AP 


P8-11-2002 








JP 


2003522956 


T 


29-07-2003 








NO 


20014972 


A 


12-12-2001 








PL 


350336 


Al 


02-12-2002 








SK 


16312001 


A3 


02-07-2002 








US 


2002191490 


Al 


19-12-2002 



US 5377104 A 27-12-1994 NONE 



US 6049508 


A 


11-04-2000 FR 


2772137 Al 


11-06-1999 






BR 


9805239 A 


03-11-1999 






CA 


2253764 Al 


08-06-1999 






DE 


69807826 Dl 


17-10-2002 






DE 


69807826 T2 


16-01-2003 






EP 


0921416 Al 


09-06-1999 






NO 


985718 A 


09-06-1999 



Form PCT/lSA/210 (patent family annex) (July 1992) 



RAPPORT DE RECHEBCHE INTERNATIONALE 



Demj 

PC 



T7mW3, 



itionale No 

'3/02874 



B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 


Documentation mintmale consultee (systeme de classification suivi des symbotes de c 

CIB 7 G01V 


dassement) 


Documentation consultee autre que la documentation minimate dans la mesure 0C1 ce 


s documents retevent des domaines sur lesquels a porte la recherche 



A. CLASSEMENT DE L'OBJET DE LA DEMANDE 
CIB 7 G01V1/28 



Selon la classification Internationale des brevets (CIB) ou a la fojs selon la classification natlonale et la CIB 



Base de donnees electronlque consultee au cours de la recherche Internationale (nom de la base de donnees, et si realisable, termes de recherche utilises) 

EPO-Internal , WPI Data, INSPEC, COMPENDEX 



C. DOCUMENTS CONSIDERES COM ME PERTINENTS 



Categorie ° Identification des documents cites, avec, te cas echeant, P'mdicatlon des passages pertinents 



no. des revendications visees 



US 6 442 489 Bl (S0K0LIN HAIM ET AL) 
27 aoQt 2002 (2002-08-27) 
colonne 3, ligne 31 - ligne 48 
colonne 4, ligne 40 - ligne 54 
colonne 5, ligne 37 -colonne 6, ligne 25 

FR 2 805 051 A (GEOPHYSIQUE CIE GLE) 
17 aoOt 2001 (2001-08-17) 
cite dans la demande 
revendi cation 1 



A 



US 5 377 104 A (SORRELLS GORDON G 
27 decembre 1994 (1994-12-27) 
abrege 



ET AL) 



US 6 049 508 A (DEFLANDRE JEAN-PIERRE) 
11 avril 2000 (2000-04-11) 
revendi cation 1 



1,7 
1,7 



| | Voir la suite du cadre C pour la fin de ta ilste des documents 



Las documents de families de brevets sont indiques en annexe 



° Categories speclales de documents cites: 

"A" document definissant Petal general de la technique, non 
considere cornme particulierement pertinent 

■E' document anterieur, mats publie a la date de depot international 
ou apres cette date 

•L' document pouvant jeter un doute sur une revendication de 
priortte ou cite pour determiner la date de publication d'une 
autre citation ou pour une ralson speciale (telle qulndiquee) 

'O' document se referant a une divulgation orale, a un usage, a 
une exposition ou tous autres moyens 

•P" document publie avant la date de depot international, mats 
posterieurement a la date de priortte revendiquee 



T document ulterieur publie apres la date de d£pot international ou la 
date de priorite et n*ap parte nenant pas a Fetal de la 
technique pertinent, mais cite pour comprendre le principe 
ou latheorie constituent la base de Pinvention 

"X" document particulierement pertinent; Plnven tlon revendiquee ne peut 
etre considered comme nouveile ou comme impliquant une activite 
inventive par rapport au document considere Isolement 

"Y" document particulierement pertinent; PInven tlon revendiquee 
ne peut etre consideree comme impliquant une activite inventive 
lorsque le document est associe k un ou plusieurs autres 
documents de meme nature, cette com bin at son etant evldente 
pour une personne du metier 

'&■ document qui fait partie de la meme famiile de brevets 



Date a laqueQe la recherche intemationale a ete effectivement achevee 



19 fevrier 2004 



Date d*expedition du present rapport de recherche Internationale 



26/02/2004 



Nom et adresse postale de Padmlnlstration chargee de la recherche Intemationale 
Office Europeen des Brevets, P.B. 5618 Patenttaan 2 
NL-2280HVRIjswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax (+31-70) 340-3016 



Fonctionnaire autorise 



Lome, B 



Formulaire PCT/1SA/210 (deuxifeme feuiDe) Quillet 1992) 



RAPPORT DE RECHERCHE INTERNATIONALE 



R^E 



Den ; ^^fcaatlonole No 

PCT7nR/02874 



Document brevet cit§ 


j Date de 


Membre(s) de la 


Datede 


au rapport de recherche 


I publication 


famllle de brevet(s) 


publication 



US 6442489 



Bl 



27-08-2002 AUCUN 



FR 2805051 A 17-08-2001 



FR 


2805051 


Al 


17-08-2001 


BR 


0104488 


A 


08-01-2002 


CA 


2366550 


Al 


16-08-2001 


CN 


1363046 


T 


07-08-2002 


CZ 


20013702 


A3 


13-02-2002 


EA 


3029 


Bl 


26-12-2002 


EP 


1173781 


Al 


23-01-2002 


WO 


0159481 


Al 


16-08-2001 


HU 


0202591 


A2 


28-11-2002 


JP 


2003522956 


T 


29-07-2003 


NO 


20014972 


A 


12-12-2001 


PL 


350336 


Al 


02-12-2002 


SK 


16312001 


A3 


02-07-2002 


US 


2002191490 


Al 


19-12-2002 



US 5377104 A 27-12-1994 AUCUN 



US 6049508 


A 


11-04-2000 FR 


2772137 Al 


11-06-1999 






BR 


9805239 A 


03-11-1999 






CA 


2253764 Al 


08-06-1999 






DE 


69807826 Dl 


17-10-2002 






DE 


69807826 T2 


16-01-2003 






EP 


0921416 Al 


09-06-1999 






NO 


985718 A 


09-06-1999 



Formula] re PCT/ISA/21 0 (annexe families de brevets) Quillet 1992) 



